Ginzburg-Landau theory of vortices in a multi-gap superconductor 
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The Ginzburg-Landau functional for a two-gap superconductor is derived within the weak-coupling 
BCS model. The two-gap Ginzburg-Landau theory is, then, applied to investigate various magnetic 
properties of MgB2 including an upturn temperature dependence of the transverse upper critical field 
and a core structure of an isolated vortex. Orientation of vortex lattice relative to crystallographic 
axes is studied for magnetic fields parallel to the c-axis. A peculiar 30°-rotation of the vortex lattice 
with increasing strength of an applied field observed by neutron scattering is attributed to the 
multi-gap nature of superconductivity in MgB2. 
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I. INTRODUCTION 

Superconductivity in MgB2 discovered a few years 
agoi has attracted a lot of interest both from funda- 
mental and technological points of view.- Unique phys- 
ical properties of MgB2 include Tc = 39 K, the high- 
est among s-wave phonon mediated superconductors, 
and the presence of two gaps Ai « 7 meV and A2 w 
2.5 meV evidenced by the scanning tunneling^ and the 
point contact^-® spectroscopies and by the heat capacity 



measurements 
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The latter property brings back the 
concept of a multi-gap superconductivity^^'^^ formulated 
more than forty years ago for materials with large dispar- 
ity of the electron-phonon interaction for different pieces 
of the Fermi surface. 

Theoretical understanding of normal and supercon- 
ducting properties of MgB2 has been advanced in the 
direction of first-principle calculations of the electronic 
band structure and the electron-phonon interaction, 
which identified two distinct groups of bands with large 
and small superconducting gapsiiii^iiSiiSiiLiSiiS Quan- 
titative analysis of various thermodynamic and trans- 
port properties in the superconducting state of MgB2 
was made in the framework of the two-band BCS 
modeh^°i^-^i^^i^^i^'^i^^i^^i^''i^^ An outside observer would 
notice, however, a certain lack of effective Ginzburg- 
Landau or London type theories applied to MgB2. This 
fact is explained of by quantitative essence of the dis- 
cussed problems, though effective theories can often give 
a simpler insight. Besides, new experiments constantly 
raise different types of questions. For example, recent 
neutron diffraction study in the mixed state of MgB2 
has found a strange 30°-reorientation of the vortex lat- 
tice with increasing strength of a magnetic field applied 
along the c-axisi^ Such a transition represents a marked 
qualitative departure from the well-known behavior of 
the Abrikosov vortex lattice in single-gap type-II super- 
conductors. Nature and origin of phase transitions in 
the vortex lattice are most straightforwardly addressed 
by the Ginzburg-Landau theory. 

In the present work we first derive the appropriate 
Ginzburg-Landau functional for a two-gap superconduc- 



tor from the microscopic BCS model. We, then, in- 
vestigate various magnetic properties of MgB2 using 
the Ginzburg-Landau theory. Our main results include 
demonstration of the upward curvature of Hc2{T) for 
transverse magnetic fields, investigation of the vortex 
core structure, and explanation of the reorientational 
transition in the vortex lattice. The paper is organized 
as follows. Section 2 describes the two-band BCS model 
and discusses the fit of experimental data on the tempera- 
ture dependence of the specific heat. Section 3 is devoted 
to derivation of the Ginzburg-Landau functional for a 
two-gap weak-coupling superconductor. In Section 4 we 
discuss various magnetic properties including the upper 
critical field and the structure of an isolated vortex. Sec- 
tion 5 considers the general problem of an orientation 
of the vortex lattice in a hexagonal superconductor in 
magnetic field applied parallel to the c-axis and, then, 
demonstrates how the multi-gap nature of superconduc- 
tivity in MgB2 determines a reorientational transition in 
the mixed state. 



II. THE TWO-BAND BCS MODEL 

A. General theory 

In this subsection we briefly summarize the thermo- 
dynamics of an s-wave two-gap superconductor with the 
aim to extract subsequently microscopic parameters of 
the model from available experimental data for MgB2. 
We write the pairing interaction as 
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(1) 

where n = 1,2 is the band index. A real space represen- 
tation is obtained from a general momentum-space 
form of the modeUii^ under assumption of weak mo- 
mentum dependence of the scattering amplitudes 
We also assume that the active band has the strongest 
pairing interaction gn = gi compared to interaction in 
the passive band 922 = .92 and to interband scattering 
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of the Cooper pairs gi2 = (721 = .93- Defining two gap 
functions 



A„(x) = -^g„„,(*„,|(x)*„/t(x)> 



(2) 



the total Hamiltonian is transformed to the mean-field 
form 



(3) 



^nix)^l^ix)^i^ix) + h.C. 



h{x) being a single-particle Hamiltonian of the normal 
metal. The constant term is a quadratic form of anoma- 
lous averages {'^ n[{x)'^ n]{x)) ■ Using Eq. (0) it can be 
expressed via the gap functions 

Econst = ^ j dx[52|Ai|2-H5-i|A2p -53(AtA2 + A;Ai)] 

(4) 

with G = det{(7„„'} = gig2 — g^- The above expres- 
sion has to be modified for G — 0. In this case the two 
equations ||2J) are linearly dependent. As a result, the ra- 
tio of the two gaps is the same for all temperatures and 
magnetic fields A2(a;)/Ai(a;) =53/51, while the constant 
term reduces to i?const — J dx jAip/^i. 

The standard Gorkov's technique can then be applied 
to derive the Green's functions and energy spectra in uni- 
form and nonuniform states with and without impurities. 
In a clean superconductor in zero magnetic field the two 
superconducting gaps are related via the self-consistent 
gap equations 
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with dimensionless coupling constants A„„' = gnn'N„i, 
Nn being the density of states at the Fermi level for 
each band. The transition temperature is given by 
Tc = (2a;£ie'^/7r)e^^/^, where ujo is the Debay frequency, 
C is the Euler constant and A is the largest eigenvalue of 
the matrix A„„/ : 



A = (An + A22)/2 + v/(Aii - A22)V4 + A12A21 . (6) 

Since A > An, the interband scattering always increases 
the superconducting transition temperature compared to 
an instability in a single-band case. The ratio of the two 
gaps at T = Tc is A2/A1 = A2i/(A - A22). At zero 
temperature the gap equations ISJ are reduced to 



A„ = ^ A„„' A„/ In . 



(7) 



By substituting A„ = 2a;u7'„e the above equation is 
transformed to 



(8) 



For 1/A ^ lnr„, one can neglect logarithms on the right 
hand side and obtain for the ratio of the two gaps the 
same equation as at T = Tc implying that A2/A1 is 
temperature independent."^'' This approximation is valid 
only for r„ ~ 1, i.e., if all the coupling constants A„„/ 
have the same order of magnitude. (For g^ — 17132 the 
above property is an exact one: the gap ratio does not 
change neither with temperature nor in magnetic field.) 
However, for 53 <C 52 < 5i, the passive gap A2 is signifi- 
cantly smaller than the active gap Ai and r2 ^ 1 so that 
the corresponding logarithm cannot be neglected. It fol- 
lows from Eq. © that the ratio A2/A1 increases between 
T = Tc and T = for small 33. Such variations become 
more pronounced in superconductors with larger values 
of A, which are away from the extreme weak-coupling 
limit A ^ 1. Ab-initio calculations indicate that MgB2 
has an intermediate strength of the electron-phonon cou- 
pling with Ai2(2i) ^ All ^ 1, making this superconduc- 
tor an ideal system to observe effects related to variations 
of the ratio of two gaps. 

The jump in the specific heat at the superconducting 
transition can be expressed analytically asiSi^Si^ 
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(9) 



where the limit T ^ Tc has to be taken for the ratio of the 
two gaps. The specific heat jump is always smaller than 
the single-band ECS result AC/C = 12/7C(3) « 1.43, 
unless Ai = A2. 



B. Fit to experimental data 

One of the striking experimental evidences of the 
double-gap behavior in MgB2 is an unusual tempera- 
ture dependence of the specific heat with a shoulder-type 
anomaly around 0.25Tcii2'2iili We use here the multi- 
band BCS theory to fit the experimental data for C(T). 
The Fermi surface in MgB2 consists of four sheets: two 
nearly cylindrical hole sheets arising from quasi two- 
dimensional px.y boron bands and two sheets from three- 
dimensional bonding and antibonding bands. ^'^''^^ The 
electronic structure of MgB2 is now well understood from 
a number of density- functional studies jiii^ji^tiSiiLiSiiS 
which generally agree with each other, though differ in 
certain details. Specifically, we choose as a reference the 
work of Kong et al.f& where the tight-binding fits for 
all Fermi surface sheets in MgB2 are provided. Using 
these fits we have calculated various Fermi surface aver- 
ages for each band. The density of states at the Fermi 
level is A^(0) — 0.41 states/eV/cell/spin, which includes 
iV^(O) = 0.16 = 0.049 0.111 states/eV/cell/spin in 
light and heavy cr-bands and iV^(O) = 0.25 = 0.124 -1- 
0.126 states/eV/cell/spin in the two 7r-bands. Note, 
that the obtained iV7r(0) is somewhat larger than the 
number 0.205 cited by Kong et al.^ while the results 
for the (T-bands agree. Because of a strong mismatch 
in the electron-phonon coupling between two group of 
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FIG. 1: Theoretical dependence of the specific heat in the 
two-band BCS model. Numbers for each curve indicate values 
of gi, 32, and ga [N-i — 0.4, N2 = 0.6). Open circles are the 
experimental data.^'^" 



I2g^jjjg^i5ji6ji7ji8 ^jjg cr-bands can be represented as a 
single active band, which has A^i = 0.4A^(0) of the total 
density of states and drives superconducting instability, 
whereas a combined 7r-band contributes N2 — 0.6A^(0) 
to the total density of states and plays a passive role in 
the superconducting instability. The above numbers are 
consistent with iVi = 0.45iV(0) and iVi = 0.42A^(0) for 
the partial density of states of the of the electrons in the 
cr-bands obtained in the other studies>i^iii 

The gap equations © have been solved self- 
consistent ly for N2/N1 = 1.5 and various values of cou- 
pling constants. The specific heat is calculated from 
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where Enk = \/ £k + is a quasiparticle energy for 
each band and np{e) is the Fermi distribution. Figure ^ 
shows two theoretical fits to the experimental data of 
Geneva group2iiS using a weak giNi = 0.4 and a moder- 
ate giNi = 0.8 strength of the couphng constant in the 
active band. Constants (72 and 53 have been varied to 
get the best fits. In the first case the gap ratio changes 
in the range A1/A2 =3.-2.5 between T = and T = 0, 
while in the second case A1/A2 ~ 2.7. Both theoretical 
curves reproduce quite well the qualitative behavior of 
C(T). Somewhat better fits can be obtained by increas- 
ing the partial density of states in the cr-band. Quantita- 
tive discrepancies between various theoretical fits and the 
experimental data are, however, less significant than dif- 
ferences between different samplesiiS We, therefore, con- 
clude that though the specific heat data clearly agree 
with the two-gap superconducting model in the regime 
of weak interband interaction, a unique identification of 
coupling constants is not possible from available data. 



III. THE GINZBURG-LANDAU FUNCTIONAL 

We use the microscopic theory formulated in the pre- 
vious section to derive the Ginzburg-Landau functional 
of a two-gap superconductor. In the vicinity of Tc the 
anomalous terms in the mean-field Hamiltonian (0 are 
treated as a perturbation Va- Then, the thermodynamic 
potential of the superconducting state is expressed as 



ri., — En 



In / T.r exp 



Va{T)dT 



(11) 



where /3 = 1/T. Expansion of Eq. Hll|l in powers of 
Va yields the Ginzburg-Landau functional. Since the 
normal-state Green's functions are diagonal in the band 
index, the Wick's decoupling of Va in Vis does not pro- 
duce any mixing terms between the gaps. As a result, the 
weak-coupling Ginzburg-Landau functional has a single 
Josephson-type interaction term: 



= dx 



ai|Ai|2 + a2|A2p-7(AtA2 



A;Ai) 



i/3l|Ai|4-f i/32|A2|4+ Ku\'^^A,\^+K2^\V,A2\ 



V. = 9. + .J^A„ ar,2 = ^-A^i,2ln^, (12) 



7C(3)jV„ 33 
^" 16^2^2 ' ^ G" 



K„ 



7C(3)jVn 
167r2T2 



$0 being the flux quantum. For 7 > 0, the interaction 
term favors the same phase for the two gaps. For 7 < 0, 
if, e.g., the Coulomb interactions dominate the interband 
scattering of the Cooper pairs and 53 < 0, the smaller gap 
acquires a 7r-shift relative to the larger gapiSiM 

The gradient term coefficients depend in a standard 
way on the averages of Fermi velocities "Vpn over vari- 
ous sheets of the Fermi surface. Numerical integration 
of the tight-binding fitsiS yields the following results: for 
the cr-band (vj,^) = 2.13 (3.55, 1.51) and 



(0.05, 0.05); for the 7r-band 



= 0.05 
1.51 (1.47, 1.55) and 
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{vj,J = 2.96 (2.81, 3.10) in units of lO^^ cm2/s2, num- 
bers in parentheses correspond to each of the constituent 
bands. The effective masses of the quasi two-dimensional 
cr-band exhibit a factor of 40 anisotropy between in- 
plane and out of plane directions. In contrast, the three- 
dimensional TT-band has a somewhat smaller mass along 
the c-axis. Using N2/N1 — 1.5 we find that the in-plane 
gradient constants for the two bands are practically the 
same K2±/Ki± « 1.06, while the c-axis constants differ 
by almost two orders of magnitude K2z/Kiz ~ 90. 

A very simple form of the two-gap weak-coupling 
Ginzburg-Landau functional is somewhat unexpected. 
On general symmetry grounds, there are possible 
various types of interaction in quartic and gradient 
terms between two superconducting condensates of the 
same symmetry, which have been considered in the 
literature The above form of the Ginzburg- 
Landau functional is, nevertheless, a straightforward 
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extension of the well-known result for unconventional 
superconductors. For example, the quartic term for 
a momentum-dependent gap is |A(k)|'' in the weak- 
coupling approximation^^i^ In the two-band model 
A(k) assumes a step-like dependence between different 
pieces of the Fermi surface, which immediately leads to 
the expression H12|l . 

The Ginzburg-Landau equations for the two-gap su- 
perconductor, which are identical to those obtained from 
Eq. H12I) . have been first derived by an expansion of 
the gap equations in powers of Ai^ Recently, a simi- 
lar calculation has been done for a dirty superconductor, 
with only intrahand impurity scattering and the corre- 
sponding form of the Ginzburg-Landau functional has 
been guessed, though with incorrect sign of the coupling 
term.-iS, Here, we have directly derived the free energy 
of the two-gap superconductor. The derivation can be 
straightforwardly generalized to obtain, e.g. higher-order 
gradient terms, which are needed to find an orientation of 
the vortex lattice relative to crystal axis (see below). We 
also note that strong-coupling effects, e.g., dependence of 
the pairing interactions on the gap amplitudes, will pro- 
duce other, generally weaker, mixing terms of the fourth 
order in A. The interband scattering by impurities can 
generate a mixing gradient term as well. 

Finally, for G = (9192 — 93) < a number of spurious 
features appears in the theory: the matrix A„„' and the 
quadratic form Q acquire negative eigenvalues, while a 
formal minimization of the Ginzburg-Landau functional 
(|12|1 leads to an unphysical solution at high tempera- 
tures. Sign of A2/A1 for such a solution is opposite to 
the sign of 53. The origin of this ill-behavior lies in the ap- 
proximation of positive integrals on the right-hand side 
of Eq. ISJ by logarithms, which can become negative. 
Therefore, negative eigenvalues of A„„/ and E'const yield 
no physical solution similar to the case when the BGS 
theory is applied to the Fermi gas with repulsion. The 
consequence for the Ginzburg-Landau theory H12I) is that 
one should keep the correct sign of A2/A1 and use the 
Ginzburg-Landau equations, i.e., look for a saddle-point 
solution rather than seeking for an absolute minimum. 



IV. THE TWO-GAP GINZBURG-LANDAU 
THEORY 

In order to discuss various properties of a two-gap su- 
perconductor in the framework of the Ginzburg-Landau 
theory we write ai = —ait with ai = Ni,t = hi{Ti/T) « 
(l-T/Ti) and Ti = (2wDe'^/7r)e~92/^^i for the first ac- 
tive band and a2 = a2o — a2t with 02 = N2, a2o = 
(All — ^22)/GNi for the passive band. 



is found from diagonalization of the quadratic form in 
Eq. Ua): 



tr 
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(13) 



For small 7, one finds tc ~ —7^/ (aia2o)- Negative sign of 
tc means that the superconducting transition takes place 
above Ti, which is an intrinsic temperature of supercon- 
ducting instability in the first band. The ratio of the two 
gaps p = A2/A1 = 7/(ck20 ~ CL2tc)- Below the transition 
temperature, the gap ratio p obeys 



a2P - 7 + ^P^(«i* + 7P) = 
Pi 



(14) 



For small 7, one can approximate p « 7/02 and due to a 
decrease of a2 with temperature, small to large gap ratio 
p increases away from tc. 



B. The upper critical field 

1. Magnetic field parallel to ttie c-axis 

Due to the rotational symmetry about the c-axis, the 
gradient terms in the a-b plane are isotropic and can 
be described by two constants Kn_L = Kn. The lin- 
earized Ginzburg-Landau equations describe a system of 
two coupled oscillators and have a solution in the form 
Ai = CQfo{x) and A2 — dofolx), where fQ{x) is a state 
on the zeroth Landau level. The upper critical field is 
given by Hc2 = /ic2'&o/27r 



hc2 = 



ait 
2KI 




K1K2 



The ratio of the two gaps p = do/co along the upper 
critical line is 

P= ■ (16) 

The above expression indicates that an applied mag- 
netic field generally tends to suppress a smaller gap. 
Whether this effect overcomes an opposite tendency to 
an increase of A2/A1 due to a decrease of a2 with 
temperature depends on the gradient term constants. 
For example, in the limit 7 <C a2o we find from l|16() 
p « 7/[q;2o — (^2 — 0,1X2/ Ki)t]. If K2 is significantly 
larger than Ki , while a2 — oi , the smaller gap is quickly 
suppressed along the upper critical field line. However, 
for MgB2 one finds K2/K1 « 1 for in-plane gradient 
terms. Therefore, the ratio A2/A1 continues to grow 
along Hc2{T), though slower than in zero field. 



A. Zero magnetic field 

For completeness, we briefly mention here the behav- 
ior in zero magnetic field. The transition temperature 



2. Transverse magnetic field 

We assume that H || y and consider a homogeneous 
superconducting state along the field direction. The gra- 
dient terms in two transverse directions x and z have 
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FIG. 2: Temperature dependence of the upper critical field 
in MgB2 for magnetic field in the basal plane. Solid line: the 
two-gap Ginzburg-Landau theory with parameters given in 
the text, squares: experimental data by Lyard et al.'^" The 
inset shows variation of the gap ratio along the Hc2{T) line 
for the same set of parameters. 



different stiffness constants Kn and Knz^ respectively. 
In a single band case, rescaling x ~* x(Kxli^z)^^^ and 
z z[Kzl i^x)^^^ allows to reduce an anisotropic prob- 
lem to the isotropic one in rescaled coordinates. A multi- 
gap superconductor has several different ratios Kn/Km 
and the above rescaling procedure does not work. In 
other words, coupled harmonic oscillators described by 
the linearized Ginzburg-Landau equations have different 
resonance frequencies. To solve this problem we follow 
a variational approach, which is known to give a good 
accuracy in similar cases. The vector potential is chosen 
in the Landau gauge A = (i7z,0,0) and we look for a 
solution in the form 



Ai 

A2 



1/4 



(17) 



where A, c, and d are variational parameters. After spa- 
tial integration and substitution X = h/ ^, h = 2ttH/^o, 
the quadratic terms in the Ginzburg-Landau functional 
become 



F2 = {-ait + hKi)\c\^ + {a2 + hK2)\4' 



(18) 



-^{c*d + d*c) 



K„ 



The determinant of the quadratic form vanishes at the 
transition into superconducting state. Transition field 
is given by the same expression as in the isotropic case 
(|15|l . where have to be replaced with Kn- The up- 
per critical field is, then, obtained from maximizing the 
corresponding expression with respect to the variational 
parameter /i. In general, maximization procedure has to 
be done numerically. Analytic expressions are possible in 
two temperature regimes. At low temperatures t 3> \tc\, 
the upper critical field is entirely determined by the ac- 



hc2 = 



ai t 



(19) 



In the vicinity of Tc, an external magnetic field has a 
small effect on the gap ratio p = d/c ~ 7/0120 and an 
effective single-gap Ginzburg-Landau theory can be ap- 
plied. The upper critical field is given by a combination 
of the gradient constants Kni weighted according to the 
gap amplitudes: 



hc2 = 



ai{t - tc) 



^iKi+p^K2)iKiz+P^K2z) 



(20) 



Since, in MgB2 one has Kiz ~ O.OIK2Z and p^ ~ 0.1, the 
slope of the upper critical field near Tc is determined by 
an effective gradient constant K^^ fa P^K2z > Kiz (while 
[Ki + p^K2) ~ Kx). Thus, the upper critical field line 
Hc2{T) shows a marked upturn curvature between the 
two regimes (|20() and (|19|l . Such a temperature behav- 
ior has been recently addressed in a number of theoret- 
ical works based on various forms of the two-band BCS 
theoryiS^iSSiSSiSi We suggest here a simpler description 
of the above effect within the two-gap Ginzburg-Landau 
theory. 

Finally, we compare the Ginzburg-Landau theory with 
the experimental data on the temperature dependence of 
the upper critical field for magnetic field parallel to the 
basal planers We choose ratios of the gradient term con- 
stants and the densities of states in accordance with the 
band structure calculations^^ and change parameters 7 
and a20, which are known less accurately, to fit the ex- 
perimental data. The best fit shown in Fig.[21is obtained 
for 0,20 1^1 — 0.65 and 7/ai = 0.4. The prominent up- 
ward curvature of Hc2{T) takes place between tc = —0.18 
(Tc = 36 K) and t ~ 0.2 (T 26 K), i.e. weU within the 
range of validity of the Ginzburg-Landau theory. The 
above values of a2o and 7 can be related to (72/51 and 
.93/51 and they appear to be closer to the second choice 
of gn used for Fig. |Q. The ratio of the two gaps, as it 
changes along the Hc2{T) line, is shown on the inset in 
Fig. El It varies from A1/A2 « 2.3 near = 36 K 
to A1/A2 « 45 at T 18 K, where the Ginzburg- 
Landau theory breaks down. Due to a large difference 
in the c-axis coherence lengths between the two bands, 
the smaller gap is quickly suppressed by transverse mag- 
netic field. Also, the strong upward curvature of Hc2{T) 
leads to temperature variations of the anisotropy ratio 
7a„ = H^2{T)/H^2iT), which changes from 7a„ = 1.7 
near Tc to 7an = 4.3 at T = 18 K. These values are again 
consistent with experimental observations}^ as well as 
with theoretical studiesi^'^1^^1^^ 



C. Structure of a single vortex 

The structure of an isolated superconducting vortex 
parallel to the c-axis has been studied in MgB2 by the 
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FIG. 3: Spatial dependencies of the gaps for various temper- 
atures with t = InTi/T ?a 1 - T/Ti and K2/K1 = 9. 



FIG. 4: Spatial dependencies of the gaps for various values of 
7 (given in units of ai) for t — 0.3 and K2/K1 — 9. 



scanning tunneling microscopy4^ Tunneling along the c- 
axis used in the experiment probes predominantly the 
three-dimensional 7r-band and the obtained spectra pro- 
vide information about a small passive gap. A large 
vortex core size of about 5 coherence lengths = 
yJ^[)/2TTH^2 reported and attributed to a fast sup- 
pression of a passive gap by magnetic field, whereas the 
c-axis upper critical field is controlled by a large gap in 
the (T-band^ The experimental observations were con- 
firmed within the two-band model using the Bogoliubov- 
de Gennes22: and the Usadel equations. We have, how- 
ever, seen in the previous Subsection that a 7r-gap in 
MgB2 is not suppressed near Hc2{T) for fields applied 
along the c-axis. To resolve this discrepancy we present 
here a systematic study of the vortex core in a two-gap su- 
perconductor in the framework of the Ginzburg-Landau 
theory. 

We investigate the structure of a single-quantum vor- 
tex oriented parallel to the hexagonal c-axis. The two 
gaps are parametrized as A„(r) = where 9 is 

an azimuthal angle and r is a distance from the vortex 
axis. Since the Ginzburg-Landau parameter for MgB2 
is quite large;^ k ~ 25, magnetic field can be neglected 
inside vortex core leading to the following system of the 
Ginzburg-Landau equations 

ani^n - 7^n' + Ml " K^{^',[ + ^j'Jr - QVn) = (21) 

for n = 1,2 {n' — 2, 1) and Q k, 1/r. Away from the cen- 
ter of a vortex, the two gaps approach their asymptotic 
amplitudes V'On 



^1 



'01 



ait -f 7p 



f3i 



'02 



'7/p - a2 



P2 



(22) 



with p obeying Eq. (|14(l . All distances are measured in 
units of a temperature-dependent coherence length de- 
rived from the upper critical field Eq. (|15|l . In order to 



solve Eq. H21|l numerically, a relaxation method has been 
used^ on a linear array of 4000 points uniformly set on 
a length of 80^ from the vortex center. An achieved ac- 
curacy is of the order of 10~^. 

The obtained results are shown in Figures 3-5, where 
amplitudes ipnif) are normalized to the asymptotic value 
of the large gap ^pQi. To quantify the size of the vortex 
core for each component we determine the distance c?„, 
where ipnif) reaches a half of its maximum value tpon- In 
the case of a single-gap superconductor such a distance 
is given within a few percent by the coherence length. In 
a two-gap superconductor the characteristic length scale 
for the large gap di remains close to ^, while d2 can 
substantially vary. Size of the vortex core is given by 
dy = max(2di, 2^2). 

Results for temperature dependence of the vortex core 
are presented in Fig. 13 The parameters a2o and 7 are 
taken the same as in the study of the upper critical field, 
while we choose K2/K1 = 9 in order to amplify effect 
for the small gap. As was discussed above, the equilib- 
rium ratio of the two gaps "002 /■0oi grows with decreas- 
ing temperature (increasing t) . Simultaneously, the small 
gap becomes less constrained with its interaction to the 
large gap and the half-amplitude distance d2 shows a no- 
ticeable growth. For K2 « Ki such a less constrained 
behavior of 4'2{f) at low temperatures does not lead to 
an increase of the core size because both gaps have sim- 
ilar intrinsic coherence lengths. 

This trend becomes more obvious if the coupling con- 
stant 7 is changed for fixed values of all other param- 
eters, see Fig. 0] For vanishing 7, the distance ^2 ap- 
proaches asymptotically an intrinsic coherence length in 
the passive band. This length scale depends on K2 
id2/di\-y=a ~ \/K2lK\ = 3), but is not directly related 
to an equilibrium value of the small gap: the small gap 
is reduced by a factor of 7 between 7 = 0.6 and 7 = 0.03, 
while the core size increases by 50% only. Therefore, the 
single-band BCS estimate ^2 = vf/{t^^2) for the char- 
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FIG. 5: Spatial dependencies of the gaps for various values of 
K2/K1 for t = 0.3. 



acteristic length scale of the small gap sometimes used 
for interpretation of experimental data*^ is not, in fact, 
applicable for a multi-gap superconductor. 

Finally, Fig. presents evolution of the vortex core 
with varying ratio K2/K1, where again Q!2o/ai = 0.65 
and 7/ai = 0.4. The apparent size of the vortex core 
dv ~ 2(^2 becomes about 5-6 coherence lengths for K2 
exceeding Ki by an order of magnitude. For K2/K1 ~ 1, 
which follows from the band structure calculations, the 
vortex core size does not change significantly compared 
to the standard single-gap case. These results gener- 
ally agree with the previous study^S though we conclude 
that unrealistically large values of K2 /Ki are required to 
explain the experiment. Different strength of impurity 
scattering in the two bands cannot explain this discrep- 
ancy either. It is argued that the 7r-band is in the dirty 
limit The coefficient K2 in Eq. (|12|l is accordingly re- 
placed by a smaller diffusion constant. The numerical 
results (Fig. Isj as well as qualitative consideration show 
that in such a case the core size for 4'2{f) can only de- 
crease. Note, that the spatial ansatz '0(r) ^ tanh(r/a) 
with a = ^ used to fit the experimental data^ should be 
applied with a = 1.8^ even for a single-gap superconduc- 
tor in the large n limit 44 



V. ORIENTATION OF VORTEX LATTICE 

Recent neutron scattering measurementsSS in MgB2 
for fields along the hexagonal c-axis have discovered a 
new interesting phase transition in the mixed state of 
this superconductor: a triangular vortex lattice rotates 
by 30° such that below the first transition field (0.5 T at 
T = 2 K) a nearest-neighbor direction is aligned perpen- 
dicular to the crystal a-axis, whereas above the second 
transition field (0.9 T) the shortest intervortex spacing is 
parallel to the a-axis We show in this section that such 



a peculiar behavior is determined by the two-gap nature 
of superconductivity in MgB2. 



A. Single-gap superconductor 

An orientation of the flux line lattice in tetragonal 
and cubic superconductors has been theoretically stud- 
ied by TakanakaJS, Recently, similar crystal field effects 
were found to be responsible for the formation of the 
square vortex lattices in the borocarbides.^^'"'^ The case 
of a single-gap hexagonal superconductor is treated by 
a straightforward generalization of the previous works. 
Symmetry arguments suggest that coupling between the 
superconducting order parameter and a hexagonal crys- 
tal lattice appears at the sixth-order gradient terms in the 
Ginzburg-Landau functional. For simplicity, we assume 
that gap anisotropy is negligible. Then, the six-order 
gradient terms derived from the BCS theory are 



Fa = 



C(7)iVo 



1 



327r6T6 
X (V.,VjVfcA*)(V;V„V„A 



{VFI VF j VFkVpiVFmVFn ) 



(23) 



The above terms can be split into isotropic part and 
anisotropic contribution, the latter being 



pan ^ _C(W 

^ 647r6T6 



(1-^) ((4.>-(4,» 



X A* 



- 15V^V2 



15V2 



A (24) 



>aA' 



In this expression x is fixed to the a-axis in the basal 
plane. (An alternative choice is the 6-axis.) If x and 
y are simultaneously rotated by angle about the c- 
axis, {Vx ± *Vj,)® acquires an extra factor e^^"^"^ . In the 
following we always make such a rotation in order to have 
X pointing between nearest-neighbor vortices. Periodic 
Abrikosov solutions with chains of vortices parallel to 
the X-axis are most easily written in the Landau gauge 
A = (-i/2/,0,0)4£ 

The higher order gradient terms Eq. (|24|l give a small 
factor (1 — T/Tc)"^ and can be treated as a per- 

turbation in the Ginzburg-Landau regime. The Landau 
levels expansion yields A(a;) = Cofo{x) + CQfQ{x) + 
where the coefficient for the admixed sixth Landau level 
is ce/co « ~{V6^./3)h'^e^"^Ke/K. When substituted into 
the quartic Ginzburg-Landau term, this expression pro- 
duces the following angular dependent part of the free 
energy: 



SFiip) 



^'^^■^'-hl2l3\co\'{\fo\'fSf6)co,{6^) 



3K 



(25) 



with |co|2 = K{K2 -/i)(l/ol')/(/3(l/ol^)). Spatial aver- 
aging of the combination of the Landau levels is done in 
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a standard way 

(l/olVo/e) ^ 
{\fo\T 



— =^cos(27rpnm)e 

^ n,m. 



co/o(a;) + cefeix) and A2(x) = do/o(a;) + defeix). Solu- 
tion of the linearized Ginzburg-Landau equations yields 
the following amplitudes for the sixth Landau levels: 



C6 



45 , ,2 15 
+ — 7rcr(n-m) - — 



(26) 



aiQ!2 - 7^ 



(27) 



where summation goes over all integer n and m and pa- 
rameters p and a describe an arbitrary vortex lattice.^® 
For a hexagonal lattice (p = 1/2, a = \/3/2), the nu- 
merical value of the lattice factor is (|/oP/o/6)/(|/o|'*) = 
-0.279. Hence, 8F{^p) ~ cos(6(^) and for 



with ai,2 — 0L\^2 + \2>K\ 2h. Subsequent calculations fol- 
low closely the single-gap case from the preceding sub- 
section. The angular dependent part of the free energy is 
obtained by substituting H27|l into the fourth-order terms: 



> 



{Ke > 0) the equilibrium angle is ip = Tr/2 (tt/G), SF{ip) = [Pi4ice+4)+P2dtide+dl)\ (l/or/o/e) • (28) 
which means that the shortest spacing between vortices 
in a triangular lattice is oriented perpendicular to the a- 
axis, while for the other sign of anisotropy the shortest 
side of a vortex triangle is along the a-axis. Thus, the 
Fermi surface anisotropy fixes uniquely the orientation of 
the flux line lattice near the upper critical field. 



The resulting expression can be greatly simplified if one 
uses (A2/Ai)2 ~ 0.1 as a small parameter. With accu- 
racy 0[(A2/Ai)'*] we can neglect the angular dependent 
part determined by the small gap. This yields in a close 
analogy with Eq. (|25|l the following anisotropy energy for 
the vortex lattice near Hc2 



B. Two-gap superconductor 

In a multiband superconductor effect of crystal 
anisotropy may vary from one sheet of the Fermi sur- 
face to another. We apply again the tight-binding 
representatiorii& to obtain a quantitative insight about 
such effects in MgB2. Explicit expressions for dispersions 
of the two hole a-bands are presented in the Appendix. 
Hexagonal anisotropy in the narrow cr-cylinders is en- 
hanced by a nonanalytic form of the hole dispersions. 
Combined anisotropy of the cr-band is {vpj} — 4.608, 



hl2f3i\co\\\fo\'f5h}Kfcos{6^) 
2 
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{a2 + K2h){a2 + ISi^a/i) 



(29) 



4y) - 



"Fx/ 

4.601, while for the 7r-band {v%^) = 1.514, 
1.776 in units of 10*^ (cm/s)®. According to 
the choice of the coordinate system^i^ the i-axis is par- 
allel to the fe-dircction and the y-axis is parallel to the 
a-direction in the boron plane. The above values might 
be not very accurate due to uncertainty of the LDA re- 
sults, however, they suggest two special qualitative fea- 
tures for MgB2. First, relative hexagonal anisotropy of 
the Fermi velocity vpnif) differs by almost two orders of 
magnitude between the two sets of bands. Second, cor- 
responding hexagonal terms have different signs in the 
two bands. In Appendix, we have shown that the sign 
difference is a robust feature of the tight-binding approx- 
imation and cannot be changed by a small change of the 
tight-binding parameters. 

We investigate equilibrium orientation of the vortex 
lattice in MgB2 within the two-gap Ginzburg-Landau 
theory. Anisotropic sixth-order gradient terms of the 
type (|24|) have to be added to the functional 1)12(1 sepa- 
rately for each of the two superconducting order param- 
eters. As was discussed in the previous paragraph the 
anisotropy constants have different signs Kqi > and 
Kq2 < and obey \Kqi\ <C \Kq2\- In the vicinity of the 
upper critical field the two gaps are expanded as Ai (x) = 



Despite the fact that we have omitted terms ~ dgde, the 
Fermi surface anisotropy of the second band still con- 
tributes to the effective anisotropy constant K^^ via lin- 
earized Ginzburg-Landau equations. Along the upper 
critical line this contribution decreases suggesting the fol- 
lowing scenario for MgB2. 

In the region near Tc the second band makes the largest 
contribution to K^^: a small factor 7^/02 ~ 0.1 is out- 
weighed by \Kqi/Kq2\ < 0.1. As a result, Kg^ is neg- 
ative and if — 0, which means that the shortest inter- 
vortex spacing is parallel to the 6-axis. At lower tem- 
peratures and higher magnetic fields the second term in 
Kq^ decreases and the Fermi surface anisotropy of the 
first band starts to determine the (positive) sign of K^^ . 
In this case, ip = tt/2 {tt/6) and the side of the vor- 
tex triangle is parallel to the a-axis. The very small 
1^61/^62! — 1-8 • 10~2, which follows from the band 
structure data^ is insufhcient to have such a reorienta- 
tion transition in the Ginzburg-Landau region. Absolute 
values of anisotropy coefficients are, however, quite sensi- 
tive to the precise values of the tight-binding parameters 
and it is reasonable to assume that experimental values of 
Ken are such that the reorientation transition is allowed. 

The derived sequence of the orientations of the flux line 
lattice in MgB2 completely agrees with the neutron scat- 
tering data^S though we have used a different scan line in 
the H-T plane in order to to demonstrate the presence of 
the 30°-orientational transformation, see Fig. 6. Condi- 
tion Kf^ = 0, or similar one applied to Eq. H28|) . defines 
a line H*{T) in the H~T plane, which has a negative 
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(/oVD/d/oH^ = 0.804. Thus, the twelveth-order har- 
monics in Eq. H3U|I has a positive coefficient and trans- 
formation between the low-field state with ip ~ Q and the 
high- field state ip = n/G goes via a phase with intermedi- 
ate values of ip separated by two second order transitions. 

The anisotropy terms of the type (|24|l also produce a 
six-fold modulation of the upper critical field in the basal 
plane. Sign of the corresponding modulations of Hc2{'p) 
should also change at a certain temperature, which is de- 
termined by a suppression of the small gap in transverse 
magnetic field and is not, therefore, related to the inter- 
section point of Hc2{T) and H*{T) lines on the phase 
diagram for -ff || c. Fig. 6. 



FIG. 6: Phase diagram of MgB2 for fields parallel to the c- 
axis. Shaded region corresponds to intermediate orientations 
of the vortex lattice separated by dotted lines of the second- 
order transitions. Dashed lines indicate the scans used in the 
experiment (vertical) and in the presented theory. 



slope at the crossing point with Hc2{T). The six- fold 
anisotropy for the vortex lattice vanishes along H*{T) 
and all orientations with different angles p become de- 
generate in the adopted approximation. The sequence of 
orientational phase transition in such a case depends on 
weaker higher-order harmonics. One can generally write 



5F{ip) ^ Ke cos(6(^) + K12 cos{l2ip) 



(30) 



where the higher-order harmonics comes with a small co- 
efficient \Ki2\ ^ l^^el- Depending on the sign of K12 
transformation between low-field (p = and high-field 
p = tt/Q {n/2) orientations, when Kq changes sign, goes 
either via two second-order transitions {K12 > 0) or via 
single first-order transition {K12 < 0). In the former 
case the transitions take place at Kq — ±47^i2, whereas 
in the latter case the first-order transition is at Kg — 0. 
These conclusions are easily obtained by comparing the 
energy of a saddle-point solution cos(6</j) = —Ke/{4:Ki2) 
for Eq. which is SFsp = -Kl/{8Ki2), to the ener- 

gies of two extreme orientations. 

In order to determine sign of the higher-order harmon- 
ics for a two-gap superconductor we expand the fourth- 
order terms in the Ginzburg-Landau functional (|I2|) to 
the next order 



6F'{p) = - [f3,cl{cl + cf) + Mlidl + df)\ if^'fl) . 

(31) 

These terms are responsible for a cos(12iy9) anisotropy 
introduced before. Similar angular dependence is also 
induced by higher-order harmonics of the Fermi veloc- 
ity vf{p), though our estimate shows that even for 
the TT-bands corresponding modulations are very small4^ 
Sign of cos(12</3) term in Eq. H31|) depends only on a 
geometric factor, spatial average of the Landau levels 
wave functions. We find for a perfect triangular lattice 



VI. CONCLUSIONS 

We have derived the Ginzburg-Landau functional of a 
two-gap superconductor within the weak-coupling BCS 
theory. The functional contains only a single interac- 
tion term between the two superconducting gaps (con- 
densates). This property allows a meaningful analysis 
of various magnetic properties of a multi-gap supercon- 
ductor in the framework of the Ginzburg-Landau theory. 
Apart from confirming the previous results on an unusual 
temperature dependence of the transverse upper critical 
field in MgB2, we have presented detailed investigation 
of the vortex core structure and have shown that the 
orientational phase transitions observed in the flux line 
lattice in MgB2 is a manifestation of the multi-band na- 
ture of superconductivity in this material. The proposed 
minimal model for the 30°-rotation of the vortex lattice 
includes only anisotropy of the Fermi surface. An addi- 
tional source of six-fold anisotropy for the vortex lattice 
can arise from angular dependence of the superconduct- 
ing gap. It was argued that the latter source of (four- 
fold) anisotropy is essential for physics of the square to 
distorted triangular lattice transition in the mixed state 
of borocarbides.^*^ For phonon-mediated superconductiv- 
ity in MgB2, the gap modulations should be quite small, 
especially for the large gap on the narrow cr-cylinders 
of the Fermi surface. Experimentally, the role of gap 
anisotropy can be judged from the position of H*{T) line 
in the H-T plane. H*(T) does not cross Hc2{T) line in 
scenarios with significant gap anisotropy^ A further in- 
sight in anisotropic properties of different Fermi surface 
sheets in MgB2 can be obtained by studying experimen- 
tally and theoretically the hexagonal anisotropy of the 
upper critical field in the basal plane. 
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APPENDIX: ANISOTROPY IN a-BANDS 

We give here expressions for the dispersions and Fermi 
surface anisotropies in the two cr-bands, which are de- 
rived from the tight-binding fits of Kong et ali^ The 
in-plane Px.y boron orbitals in MgB2 undergo an sp^- 
hybridization with s-orbitals and form three bonding 
bands. At kj_ = these bands are spht into a non- 
degenerate A-symmetric band and doubly-degenerate E- 
symmetric band, which lies slightly above the Fermi level. 
Away from the k_L — 0-line the i?-band splits into light 
and heavy hole bands. Their dispersions obtained by ex- 
pansion of the tight-binding matrijii^ in small k± are 



ez(k) = e(fc^)-2ti 



1 1-2 



dk 



4 2ff(k) + 1 - 1/d 

384(1 + d) 
4 2g{k) + 7 + 9d 



8 -L 384(1 -Hd) 



where e{kz) = Sq — 2tz cos k^ and ^(k) = (fc^ — I5k^ky + 
15k^ky — ky)/k^,. The tight-binding parameters pre- 
sented in Rcf. Ig are Eq = 0.58 eV, t± = 5.69 eV, 
tz = 0.094 eV, and d = 0.16. The six-fold anisotropy is 
given by unusual nonanalytic terms, which are formally 
of the fourth order in k. Appearance of such nonanalytic 
terms is a direct consequence of the degeneracy of the two 
bands at fc = 0. For example, a nonanalytic form of e(k) 
is known for four-fold degenerate hole bands of Si and 
Gcj^ which have cubic anisotropy already in 0{k^) or- 
der. Nonanalyticity of e;,/i(k) leads to a relative enhance- 
ment of the hexagonal anisotropy on two narrow Fermi 
surface cylinders. This anisotropy has opposite sign in 
light- and heavy-hole bands. The net anisotropy of the 
combined u-band is determined mostly by the light-holes, 
which have larger in-plane Fermi velocities. 
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